Fermi surface nesting and the origin of Charge Density Waves in metals 
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The concept of a CDW induced by Fermi-surface nesting originated from tiie Peierls idea of 
electronic instabilities in purely ID metals and is now often applied to charge ordering in real low- 
dimensional materials. The idea is that if Fermi surface contours coincide when shifted along the 
observed CDW wave vector, then the CDW is considered to be nesting-derived. We show that in 
most cases this procedure has no predictive power, since Fermi surfaces either do not nest at the right 
wave vector, or nest more strongly at the wrong vector. We argue that only a tiny fraction, if any, of 
the observed charge ordering phase transitions are true analogues of the Peierls instability because 
electronic instabilities are easily destroyed by even small deviations from perfect nesting conditions. 
Using prototypical CDW materials NbSe2, TaSe2, and CeTes, we show that such conditions are 
hardly ever fulfilled, and that the CDW phases are actually structural phase transitions, driven by 
the concerted action of electronic and ionic subsystems, i.e., q-dependent electron-phonon coupling 
plays an indispensable part. We also show mathematically that the original Peierls construction 
is so fragile as to be unlikely to apply to real materials. We argue that no meaningful distinction 
between a CDW and an incommensurate lattice transition exists. 
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One common misconception in modern solid state 
physics is that Fermi surface (FS) nesting is always, 
or nearly always responsible for a charge density wave 
(CDW) . While many materials experience a structural or 
magnetic transition with a wave vector that is incommen- 
surate or badly commensurate with the high-symmetry 
phase, and while in some cases a visual inspection of the 
FS seems to reveal nesting parts with roughly the same 
wave vector, a quantitative search aimed at finding a 
nesting-driven instability at the experimental CDW vec- 
tor on the level of one-electron energies practically always 
iailsLL^dn^. This failure is a symptom of a larger mis- 
conception about CDW's, specifically, that they are the 
result of a purely electronic instability along the lines of 
the Peierls instability in one dimension (ID), and an even 
bigger misconception that the structure of the electronic 
susceptibility in the reciprocal space can be revealed by 
inspecting the Fermi surface alone, without analyzing the 
high-energy electronic excitations. 

In the Peierls picture, lattice distortion is a sec- 
ondary effect that arises in response to an electronically 
driven charge redistribution that would occur regardless 
of whether or not the ions subsequently shift from their 
high symmetry positions. In real materials, the electronic 
and ionic instabilities always occur simultaneously. Com- 
putational attempts to stabilize a CDW without allowing 
the ions to move have failed in all cases we are aware of, 
most particularly for protototypical CDW metal NbSe2. 
We will show that the concurrence of the two transitions 
is not a coincidence but arises from the fact that CDW 
formation relies on the lattice distortion as an essential 
element and not the reverse. The necessity of strong q- 
dependent electron-phonon coupling indicates that Fermi 
surface nesting, a purely electronic effect, may help CDW 
formation, but cannot be the only driving force behind 
the CDW phenomenon^. 

In the literature the term CDW is used in two differ- 
ent senses. In some cases, a structural transition with an 
incommensurate or long period is termed a CDW, regard- 



less of its origin, while in other cases, the words CDW are 
reserved for Peierls-like instabilities that occur due to a 
divergency in the real part of the electronic susceptibility, 
so that the electronic subsystem would be unstable per 
se, even if the ions were clamped at their high symmetry 
positions. We choose the latter definition for two reasons: 
first, because many accepted CDW materials actually ex- 
hibit a commensurate CDW phased, and second, because 
it allows for a distinction between a general incommensu- 
rate lattice transition (ILT) and a CDW. Moreover, the 
archetypal CDW, the ID Peierls transition in a half-filled 
band, is commensurate with just a doubled unit cell. To 
fit the definition of a Peierls system (and therefore for 
a CDW, for the purposes of this paper) a system must 
satisfy several requirements: (a) there must be substan- 
tial nesting of the FS. Note that a quantitative measure 
of the FS nesting, sometimes called the "nesting func- 
tion" is nothing but the low-frequency limit of the imag- 
inary part of the bare electronic susceptibility, Xo(q), in 
the constant matrix element approximation^. This must, 
correspondingly, peak at the CDW wave vector, (b) the 
nesting-derived peak must carry over into the real part 
of the susceptibility, Xo(q); the same wave vector, be- 
cause it is x'(q) that defines the stability of the elec- 
tronic subsystem, (c) the peak in Xo(*l) i^ust translate 
into a divergence in the full electronic susceptibility to 
cause the electronic subsystem to be unstable even with- 
out any ionic shifts, and (d) all phonons must soften at 
the CDW wave vector, not only the one corresponding to 
the mode that eventually gives rise to the observed CDW 
(except maybe a few that cannot couple to this electronic 
instability for a particular symmetry reason). With re- 
spect to the final point, we are not aware of any material 
in which such softening for multiple modes has actually 
been observed, although in theory it is unavoidablo^. In 
fact, we intend to show that this definition as a whole 
is not fulfilled in real systems and therefore, aside from 
convention, there is nothing to distinguish a CDW from 
a structural phase transition. 
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The issue is even more confounded by the fact that 
the words "nesting" and "Peierls transition" are also in- 
terpreted differently by different researchers. We use the 
first term in its literal sense: Fermi surface parts are 
nested if, when shifted rigidly, they coincide with other 
Fermi surface parts. Some authors, appreciating the fact 
that this is not sufficient for a peak in x'(q), distinguish 
between "real nesting" that exists not only at zero fre- 
quency, but also within some finite range of transition 
energies (as we discuss later in connection with CeTes, 
this imposes an additional constraint on the Fermi ve- 
locities), and "false nesting" that has infavorable Fermi 
velocities. As for the second term, there is a tendency to 
use Peierls transition as a synonym for dimerization {e.g., 
dimerization in VO2 is often called a Peierls transitior^i^) . 
Again, wc adhere to a more strict definition that reserves 
this nomenclature for a transition driven by a lowering of 
one-electron band energy caused by the opening of a gap, 
as in the original Peierls model. Note that the doubling 
of the unit cell may occur via dimerization, but can also 
easily occur via zigzag-ing of the atomic chain or via an 
even more complicated pattern. 

In Section I, we consider the classical ID Peierls tran- 
sition, for which the conditions listed above are fully sat- 
isfied and in which FS nesting is indeed expected to give 
rise to CDW formation (the results for a nested 2D sys- 
tem, such as the nearest neighbor one-band TB model, 
are qualitatively the same). We will demonstrate how 
fragile this construction is even in ID, since a CDW of 
fully electronic origin is exponentially weak. Further- 
more, we will show how rather small deviations (on the 
order of what is expected in real materials) from the per- 
fect model are sufficiently strong to suppress FS nesting- 
driven CDW formation. In Section II, we perform den- 
sity functional calculations of several real materials that 
are commonly considered to be CDW systems and show 
that, even in these canonical systems, electron-phonon 
coupling and not nesting is at the heart of the CDW phe- 
nomenon. For clear understanding in both following sec- 
tions, it is necessary to point out that an electronic CDW 
instability is not induced by a divergence in the imagi- 
nary part of susceptibility, x"(q): which is the function 
that reflects the FS topology and can be easily measured 
experimentally by neutron scattering. Instead, it is the 
real part, x'(q), which must diverge in order to trigger 
an electronic CDW. Unfortunately, x'l£t)i is difficult to 
map experimentally (see, however, RefUll). We write the 
two parts, in the constant matrix element approximation 
as: 

/M^/(£i^ (1) 

lim x"(q,t^)/t^ = V(5(£fc - eF)<5(£fe+g - e_F)(2) 
oj— *o ^ — ^ 

k 

Since x"(qi 0) is easier to calculate, it is often pre- 
sented in first-principle studies as a quantitative test of 



the FS nesting (which it is) and/or as a gauge of a ten- 
dency to CDW formation (which it is not)i^ii^iiiii^ii^iil. 



I. A MODEL PEIERLS SYSTEM 

In the ID Peierls system, we have a parabolic band of 
noninteracting particles in a periodic external potential 
that disperses as E = (in Ry units) and has a particu- 
lar Fermi vector ^ = kp, where fi is the chemical poten- 
tial. One assumes that the band is half filled so that the 
first reciprocal lattice vector G = 2n/a = 4k p. Peierls 
was the first to point out that such a system is formally 
unstable against any doubling of the unit cell, because 
it creates an additional potential with a non-zero com- 
ponent at q = 2kF, V{q) ^ and opens a gap. He em- 
phasized that a gain in the one-electron energy depends 
on the amplitude of lattice distortion logarithmically (as 
M^logw), while the elastic energy is normally quadratic, 
and therefore the ground state would always correspond 
to a nonzero distortion. It appears that actual numer- 
ical calculations do not necessarily find a distortions^. 
Ashkenazi ei afi^ argue that when the electronic suscep- 
tibility is nonanalytic, the elastic energy may be as well 
and therefore it is not guaranteed that for an infinites- 
imal distortion, the one-electron energy will be larger. 
This same point can be argued as an inability to cleanly 
partition the total energy into the one-electron energy 
and the elastic energy in a real system, where the par- 
ticles in question are interacting electrons. Instead, the 
full expression for the total energy (in density functional 
theory, for instance) must be analyzed. Following Rcf. 

we can write, using implicit matrix notation (in real 
or reciprocal space) the change of energy arising from the 
nuclear displacement in terms of the change of the poten- 
tial of the nuclei, Wext, and the induced density change, 
6n, as: 



bEtot = -(l/2)5nx-'(5n (3) 

8n = x^Vext (4) 

X = Xo/e (5) 

e = 1 - fiXo = (1 + WiX)~^ (6) 

5Etot = -{ll2)5VextX5Vext (7) 



= -{\/2)5VextX0{l-V,XG)-HVexU (8) 

where Vi is the total DFT interaction, including the 
Coulomb and exchange-correlation kernel. Note that we 
use a sign convention such that in a stable system xo < 0- 
Neglecting the matrix character of these equations, we 
immediately observe that a divergence in xoj as found 
by Peierls, is by itself insufficient to cause an instabil- 
ity, since the total susceptibility is bounded by —Vi and 
does not diverge. This is not necessarily true if the ma- 
trix character (or local fields) is taken into account and 
such an inclusion results in a complicated formula for 
the phonon frequencies, known as the Pick-Cohen-Martin 
formula as discussed in Ref. [Toi . However, without these 
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considerations, this result demonstrates that there is no 
direct relation between the Peierls instability in a system 
of noninteracting particles and CDWs in real systems. 
The same point can be made using a linear response 
approaches: an interacting half-filled electronic system 
is stable against infinitesimal perturbation. Only a fi- 
nite distortion, and only if e-ph coupling is larger than a 
critical value, can be stable. 

This is already a very serious reservation, but nonethe- 
less it is instructive to step back and investigate the 
"classical" Peierls instability in a noninteracting system. 
There is no question that the susceptibility of this sys- 
tem is logarithmically divergent, but there are interest- 
ing, and important questions left to ask. First, is this 
divergency robust with respect to small deviations from 
a "perfect nesting" , as is always the case in real mate- 
rials? Second, where is the energy gain associated with 
this instability collected? That is, can electrons away 
from the Fermi surface be effectively neglected or must 
the effects of lower filled states be taken into account? 

The standard expression for the real part of suscepti- 
bility reads 
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FIG. 1: (color online) A comparison of x'il) under ideal ID 
conditions with perfect nesting at T=0 to x'il) under vari- 
ous non-ideal conditions. Even moderate deviations from the 
ideal (such as those found in real materials) rapidly reduce 
the divergence to a relatively weak enhancement. 
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This expression is normalized to l/kp at q (the 
overall scale is not important at the moment) and has 
a very weak, logarithmic divergency at q — ±2fci?. To 
illustrate just how weak it is, we assume a relaxation rate 
7, corresponding to the Drude relaxation rate in optics, 
and recalculate x'{q)- The new result reads: 
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(10) 



The divergency has been reduced to merely an enhance- 
ment of x'(±2fcF) over x'(0) by a factor of ln(l-h^)/4. 
For typical 7 of the order of 0.1-0.2 eV this enhancement 
is by a factor of 2-2.5. One can also add that at any finite 
temperature, even without relaxation, x'(±2fci?)/x'(0) — 

ln(2^^^j^^)/2 sa ln(|^)/2. For typical Fermi energies and 
T=:10K (most observed ILTs occur at higher tempera- 
tures) the enhancement is again only a factor of the order 
of four. 

Thus, carrier scattering and Fermi function broadening 
are suffucient to reduce the nesting-induced divergency 
to a minor structure in x'iq)- But an even more severe 
effect is caused by geometrical deviations of the Fermi 
surface from perfect nesting. A common procedure in 
the search for nesting vectors in a particular fermiology 
is to copy a quasi-2D Fermi surface cut onto transparent 
paper and slide it until some piece of the displaced Fermi 
surface visually coincides with another piece of the orig- 
inal plot. It is instructive to give a quantitative gauge of 
what constitutes a "good nesting" vs. a "bad nesting". 
Assuming that the "nested" parts really nest only up to 
some Sk in the reciprocal space, we observe, by averaging 
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which gives: 
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In the small 5k limit, x'(2fcF)/x'(0) ~ 1/2[(1 + 
hi{AkF / Sk)]. In other words, if 2D Fermi lines coincide 
within 5% of the Fermi vector, the corresponding en- 
hancement of the susceptibility is about a factor of three. 
It is also important to remember that real materials are 
quasi-2D, not 2D. Any dispersion in the third direc- 
tion of the order of Sk brings us again to EqlTT] In 
a one-band case one can estimate the ratio kp/Sk as 
~ P±/p\\- This shows that anything with a cal- 
culated transport anisotropy of less than one order of 
magnitude is 3D from a "nesting point of view" . Fig[l]il- 
lustrates the effect of various deviations from the perfect 
Peierls picture on the divergence of x'iq)- 

To further emphasize the non-role that FS nesting 
plays, let us investigate a system in which an ILT with 
Q = 2kF has actually occurred, and as a result a gap 
equal to 2V has opened, with the corresponding energy 
gain in the one-electron energy. The common wisdom 
is that energy is gained predominantly near the former 
Fermi energy, where the lowering of one-electron states 
is the largest. But, is this really so? In first order per- 
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turbation theory, the new one-electron spectrum is 



Ek = 7^ 1 (l^j 



SEk 



Ek - El 



AEk + V^Ek' + W 



(13) 



where 5Ek is the energy shift of the state \k) and AiJ^ is 
the energy difference between the two states connected 
by the nesting vector Q. An inspection of this expression 
shows that integrating it over lS.Ek docs not diverge at 
Ai?fe = (that is, at the Fermi energy), but would di- 
verge at large energies, /S.Ek oo, if the integration were 
not limited by the bandwidth. Substituting Ek = k''^, we 
get 

SEk = y'4fc|(fc - kp)^ + V^ + 2kF{k - kp). (14) 
The total energy gain, Eq = /-fe SEkdk is: 
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The first term corresponds to integrating over the region 
AEk < V, and neglecting AEk, and the second corre- 
sponds to integrating over the region AEk > V, and ne- 
glecting terms smaller than V^. In other words, the first 
(non-divergent) term comes from states near the Fermi 
level, where the gap opens up, and the second term (di- 
vergent at large /x,) comes from the rest of the states 
below the Fermi level and down to the bottom of the 
band. Since an actual instability is always a competition 
between an electronic (in this picture, one-electron) en- 
ergy gain and an elastic energy loss, it is apparently more 
important to optimize the energy gain from all occupied 
states than to open a gap over the largest possible part 
of the Fermi surface. 

Since, as we have just found, all occupied states arc 
important for a ILT, it is important to keep in mind the 
multiband character of real solids, and, correspondingly, 
the contribution of the interband transitions to x'(q). 
To get a feeling for the interband effects, let us consider 
an insulating system which, instead of the free-electron 
band discussed above, has a fully occupied first neighbor 
tight-binding band, Ei{k) = — cos(fca) separated from a 
similar empty band by an energy e, E2{k) = e + cos{ka) 
Let us assume for simplicity a constant interband dipole 
matrix element equal to unity. The susceptibility is given 
by 
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(16) 



aJe2 -4sin^(qa/2) 



This susceptibility is enhanced at q = n/a (at the edge of 
the BZ) compared to its minimum value (at g = 0) by a 
factor of e/ve^— 4. This enhancement can be very large 
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FIG. 2: (color online) A schematic showing the original en- 
ergy bands of our model system in comparison to the bands 
after a transition has caused a gap at the Fermi energy. The 
energy associated with gap formation is shaded (red) and con- 
tinues down to the bottom of the band. Representative energy 
states connected by the nesting vector are shown (AEk) along 
with infinitesimal energy gain (SEk) associated with the gap 
opening at each k. 



even for bands with a large relative shift. Using for exam- 
ple e = 2.25, (giving a gap of 1/8 of the bandwidth), the 
enhancement is 2.2, quite comparable to what one might 
expect in a realistic nesting scenario. This demonstrates 
again that the Fermi surface topology is unlikely to be a 
factor in CDW formation. Note that the band structure 
that we have used is not special in any way, as opposed 
to the half-filled band exhibiting a Peierls instability. 



II. FIRST PRINCIPLES CALCULATIONS OF 
REAL MATERIALS 

In the following sections, we present three examples of 
real materials, commonly thought to be canonical exam- 
ples of nesting-driven CDW systems: NbSe2, TaSe2, and 
CeTea. Using first principles density functional theory 
calculations, we show that FS nesting is not the driving 
force behind the observed CDW in any of the three com- 
pounds. We additionally investigate a chain of Na ions 
and show that even for this artificially perfect ID system, 
the strong FS nesting fails to produce any appreciable 
CDW when the ionic positions are fixed. This is per- 
fectly in line with our contention that electron-phonon 
coupling is necessary to instigate the ILT and that even 
in the Peierls formulation, the resulting CDW is expo- 
nentially small. 

For all investigated compounds, the first principles cal- 
culations were performed using the well-known Wien2K 
packaged! with the local density approximation (LDA) 



to the exchange correlation potential^^ . For NbSe2 and 
TaSe2, it was found that spin-orbit couphng has a finite 
effect on the band structure and Fermi surface and there- 
fore was taken into account. To get a good energy mesh, 
we calculated eigenvalues at nearly 15,000 k points in 
the full BZ. For CeTes, the partially filled /-states pose 
a well-known problem for the mean-field LDA method- 
ology by partially filling each of the / states rather than 
completely filling some and emptying all others. We have 
addressed this shortcoming by using LDA-I-U in the fully 
localized limit scheme^. We applied a U of 4.5 eV and 
successfully reproduced the state observed by magne- 
tization measurements^ for the Ce ion, simultaneously 
removing all / bands from the Fermi energy. None of the 
quantities we calculate are expected to be sensitive to the 
precise value of U and wc did not investigate the effects 
of changing it. We used a mesh of approximately 30,000 
k points in the full BZ to calculate the energy eigenvalues 
used for the susceptibilities. For all compounds, we used 
a temperature smearing of 2 mRy during numerical inte- 
gration. For calculation of the chain of Na atoms we used 
the Vienna Ab-initio Simulation (VASP) package with 
the projector augmented wave (PAW) basis se t^^i^^'^^ . 
We used several different pscudopotentials with varying 
levels of hardness for the Na core, all giving identical 
results to within the accuracy of our calculations. 



A. NbSe2 and TaSe2 

Quasi-two-dimensional NbSe2 belongs to a family of 
layered dichalcogenides that undergoes a CDW transition 
sometimes thought to be related to FS nesting. A calcu- 
lation of both the real and imaginary parts of the one- 
electron susceptibility^^ shows this to be unequivocally 
not the case. Although FS nesting does exist and does 
produce a peak in x"(q)j it is located at q = (1/3,1/3,0) 
and not at the observed qcDW = (1/3,0,0) (See Ref [l9l 
for pictures). On the other hand, a weak peak at c^cdw 
appears in x'(q). This peak is not strong enough alone 
to stimulate the observed CDW transition, but must be 
assisted by elcctron-phonon coupling at the same wave 
vector. To verify this, we performed two first principles 
calculations of NbSe2 in a supercell corresponding to the 
observed CDW vector. In one case, we clamp the ions to 
their high symmetry positions and in the other we allow 
them to shift. When allowed, the ions do indeed shift 
and an instability at the right wave vector, (1/3,0,0), is 
reproduced. In the case where the ions are forced to be 
stationary, both ionic and electronic systems remain sta- 
ble in their high symmetry states. Even if wc prompt the 
electronic system by artificially redistributing the charge 
along the known CDW direction, it relaxes back to its 
high symmetry state. Our calculations on this system 
provide a clear example of the lack of infiuence of Fermi 
surface nesting with respect to CDWs: a strong FS nest- 
ing that produces a sharp peak fails to gives rise to a 
CDW at the associated wave vector, while a peak wholly 
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FIG. 3: (color online) The Fermi surface of TaSe2 shown with 
a) no shift of the Fermi level b) a downward shift of the 
Fermi level (-40 meV). The topology of the Fermi surface, 
especially the light (yellow) sheets changes appreciably with 
a small change in Fermi energy, and brings the surface into 
good agreement with experiment 



unrelated to FS nesting but at the correct vector appears 
in the real part of the susceptibility and, in conjunction 
with electron-lattice effects (but not without!) does pro- 
duce a CDW. 

We have also investigated an isostructural material, 
TaSe2, that has an observed CDW at approximately the 
same wave vector as NbSe2. The Fermi surface of the 
Ta compound is different from that of NbSe2, as was 
most recently observed by Rossnagel et al^ who also 
point out that the measured surface differs from the cal- 
culated onoiii^i^^. This is somewhat surprising because 
calculated and observed Fermi surfaces arc in quite good 
agreement for NbSe2. One reason for the discrepancy 
between NbSe2 and TaSe2 and between theory and ex- 
periment is spin-orbit coupling, which is strong for the 
heavy Ta ion. This coupling has not been included in 
most previously published TaSe2 band structures and 
Fermi surfaces^Si^i^, with the result that they differ less 
than they should from NbSe2 band structures. The spin- 
orbit interaction non-trivially changes the band disper- 
sion near the Fermi energy, particularly along F — K and 
leads to a different Fermi surface topology. Scalar rcl- 
ativistic effects are responsible for the lower Se band in 
TaSe2 (compared to NbSe2) which removes it from the 
Fermi surface entirely. Although we include both scalar 
relativistic effect and spin-orbit coupling in our calcu- 
lations and find qualitative differences between the Nb 
and Ta di-selenides, we still find that the TaSe2 Fermi 
surface differs somewhat from the ARPES observed sur- 
face. However, nearly perfect agreement can be achieved 
by a small shift of the Fermi energy (about 0.04 eV). 
In Fig. [21 the Fermi surface both with and without the 
small shift in Fermi energy is shown. 

The susceptibilities presented here are calculated using 
the Fermi energy shift necessary to bring theoretical and 
experimental Fermi surfaces into agreement, but in fact 
this shift has no visible effect whatsoever on either the 
real or imaginary parts, despite its dramatic effect on the 
FS topology. In Fig[4l the real and imaginary parts of 
the susceptibilities for TaSe2 are presented. Looking at 
the susceptibilities in FiglH one can see some similarity 
between the peak structures in this compound and in 
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FIG. 4: (color online) The imaginary (left) and real (right) 
parts of the susceptibility for TaSe2 . The nesting peaks (imag- 
inary part) do not correspond to the observed CDW wave 
vector, while the very weak peaks in the real part do. 



NbSe2i^, most particularly in the location of the nesting- 
driven peak in x"(q) vs. the peak in x'(q)- Again we 
find that the peaks are in different locations, indicating 
that nesting cannot give rise to a charge instability. 

A recent study of these two materials^^ concluded that 
both materials have CDWs driven by strong electron- 
phonon coupling in the presence of weak nesting. Our 
calculations show that, indeed, electron-phonon coupling 
is behind the transition, but we can eliminate Fermi sur- 
face nesting from the phenomenon entirely. There is no 
nesting at all, not even weak, at the CDW wave vector. 
The electronic susceptibility structure, which is indeed 
favorable to a CDW at the right wave vector, is due to 
finite energy electronic transitions and not to a Fermi 
surface geometry. 



B. CeTea 

CeTes is another layered material, (see FigO also 
belonging to a family of compounds RTea, R = rare 
earth, all of which exhibit GDWiM^^^. The two- 
dimensionality is considerably stronger in this series than 
in the dichalcogenides discussed above, producing an eas- 
ily visible Fermi surface nesting between strongly two- 
dimensional Fermi sheets. As in the former series, the 
strongest nesting peaks seen in x"(q) produced by this 
nesting are far away from qcDW and do not carry over 
into x'(q)- However, much weaker peaks associated with 
a different nesting do appear in x"{q) at the observed 
CDW wave vector and are then strongly enhanced by 
contributions away from the Fermi energy to eventually 
produce peaks in x'(q). It is an interesting and instruc- 
tive exercise to track down the origin of both kinds of 
nesting peaks to examine why the strongest ones are ir- 
relevant and why the weaker peaks are present through- 
out the energy spectrum, finally resulting in a peak in 
x'(q). To begin, we describe how the structure of the 
CeTca very simply gives rise to its Fermi surface. CeTea 
is composed of two different types of layers: those con- 
taining staggered Ce and Te ions and those containing Te 
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FIG. 5: (color online) a) The structure of CeTes showing the 
Te (large spheres) planes interposed with Ce (small spheres) 
and Te staggered units. One of the two non-equivalent Te 
ions in the Te-only planes is shaded lighter to distinguish it. 
b) A schematic of the quasi-lD tight-binding model employed 
to illustrate the origin of the Fermi surface, c The quasi-lD 
Fermi sheets resulting from the nearest-neighbor tight binding 
model. The new BZ corresponding to the larger unit cell in 
which the two different Te ions are distinguished is shown as 
a darker (blue) diamond. 



ions only (see FiglS]). The states near the Fermi energy 
have predominantly Te p character and come from the 
pure Te layers. Assuming a nearest neighbor only tight 
binding model of Te Px and Py orbitals (See Fig[5}3) in 
these layers with w bt^^ as in FiglSb, we produce the 
crosshatched pattern of slightly warped one-dimensional 
Fermi sheets shown in Fig[St. This model and the values 
of the tight-binding parameters are very similar to those 
developed in Ref. |3^. 

Nearest neighbor Te ions are inequivalent due to sym- 
metry breaking imposed by the stacking of layers along 
the third direction. This results in a larger unit cell and a 
smaller, rotated BZ, shown superimposed in Fig|SJ:. The 
full, calculated FS is shown in Fig[6l A very good fac- 
simile of it is achieved by folding the quasi-lD sheets of 
our TB model down into the new zone. Though small 
gaps appear during the folding down process, the long 
ID ribbons are still clearly visible and both the nesting 
peak and the peak in x'il) can be traced back to these 
sheets. 

As can be seen by the cartoons in FiglH there is an ex- 
cellent (ideally perfect in the nearest-neighbor TB model) 
FS nesting between two of the quasi-lD ribbons along 
the (110) direction of the BZ. This produces very sharp 
nesting-derived peaks in x'(ci)i away from the observed 
nesting vector. A moderate peak at q = cicDW is visible, 
but by far the strongest peaks are located elsewhere in 
the Brillouin zone. Any CDW directly stemming from 
FS nesting would occur first at wave vectors correspond- 
ing to these alternate spots rather than at the observed 
cicDW- In the real part of the susceptibility, the strongest 
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FIG. 6: (color online) a) The full, calculated Fermi surface of 
CeTea (thin solid lines, green and brown online), with the one- 
dimensional bands obtained from a nearest neighbor tight- 
binding model overlayed on top (thick dashed crimson lines). 
The Fermi surface in essence can still be thought of as inter- 
secting ID Fermi sheets, even after the bands have been folded 
down into the lower symmetry cell (see text), b) A shift of 
one of the quasi-lD Fermi sheets (dashed lines represent the 
shifted FS) along the (1,1,0) direction produces nearly per- 
fect nesting, c) A shift along the observed CDW wave vector 
direction (1,0,0) produces imperfect nesting. 



FS nesting peaks are suppressed dramatically, leaving 
only a peak at the correct wave vector (Fig[6|. This is 
entirely due to an effect of the finite-energy transitions 
(note the non-zero value of terms in Eq. [T]even for widely 
spaced eigenvalues) that appear in the real part of suscep- 
tibility but not in the nesting function (imaginary part). 
The reason that the bulk of the x'(q) peak height at q 
= QC-DW comes from contributions away from the Fermi 
energy can be understood by an examination of the band 
structure: the velocities of the electronic states connected 
by HcDW (cfc and efc-i-<j) are nearly equal and opposite at 
the Fermi energy (note that for the wave vector cinest hi 
FiglSjD, the opposite is true, Vk and Vk+q have the same 
sign). Therefore nearby states -I- v6k and et+q — v5k 
are also of equal energy and are connected by (\cdw- 
These are located above and below the Fermi energy re- 
spectively and do not contribute to the nesting function 
but do contribute to the real part of the susceptibility. 
This phenomenon has been pointed out earlier in Ref. |j 
under the title of 'hidden nesting'. Obviously, if there 
are regions of the energy spectrum with equal and oppo- 
site velocities, they will contribute heavily to x'(q) even 
if they are imperfectly nested at the Fermi energy itself. 
On the other hand, strong nesting between states exactly 
at the Fermi energy may die away quickly in other parts 
of the energy spectrum and contribute little to x'(q). 
This is the essence of why well nested Fermi surfaces fail 
to produce a peak in x'(q) not only for CeTea but also 
for the dichalcogenides discussed in the previous section. 
The Fermi surface is simply too small a part of the en- 
ergy range from which x'(q) collects to have a decisive 
effect. 




FIG. 7: (color online) A diagram showing the imaginary (top) 
and real (bottom) parts of the susceptibility as a function of 
qx,qy,qs ~ 0. The arrow connects the strongest peak in the 
imaginary part to its corresponding position in the real part. 



C. Na chain 

To best approximate a perfect Peierls system, we cal- 
culated the ground state of a chain of Na ions separated 
by 38 A in the two non-chain directions. The Na-Na dis- 
tance was relaxed along the chain to its optimal value of 
3.34 A. As we saw previously with NbSe2, no CDW could 
be stabilized if the ions were clamped to their high sym- 
metry positions, even though this system is ideal in every 
respect, i.e. the electronic-only CDW is unstable even 
when no mitigating factors such as higher dimensional- 
ity or imperfect nesting are present. Surprisingly, a re- 
laxation of the ionic positions along the one-dimensional 
chain also failed to produce any deviation from the high 
symmetry, equally spaced arrangment, even when an ini- 
tial dimerization was imposed. However, when the ions 
were allowed a larger dimensional freedom, the system 
distorted into a zigzag configuration as shown in Fig. [S] 
The ions in the zigzag arrangement are 3.43 A apart and 
the angle formed by the distortion is 152 degrees. Al- 
though this distortion doubles the unit cell and lowers the 
total energy, it does not create the expected gap at Ef. 
The coupling between matrix elements at the BZ bound- 
ary (k::=±7r/a) that would give rise to a gapped system 
does not occur because the distortion is two dimensional, 
while the wavefunctions are one-diemnsional. Thus, in- 
tegration over the direction perpendicular to the chain 
in the dimension in which the distorion occurs drives the 
matrix element to zero. Precisely the same result can be 
obtained using a tight-binding formulation. Regardless of 
the number of neighbors included in the model, the states 
at the edge of the BZ remain degenerate, i.e. there is no 
gap. On the other hand, an enforced dimerization of the 
ions does produce a gap at E^^, but is unstable. The true 
origin of the two-dimensional distortion is not yet entirely 
clear, but the Peierls mechanism can be ruled out. Thus, 
under perfect conditions, no electronic-only CDW forms, 
and no electron-ion interaction assisted CDW forms ei- 
ther. This indicates that the Peierls formulation is even 
weaker than we originally set out to prove. Not only do 
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FIG. 8: The relaxed configuration of an initially one- 
dimensional chain of Na ions. If the dimensionality of the 
chain is restricted, no distortion at all occurs. 

small deviations from the ideal system destroy the diver- 
gence that is purported to cause CDW formation, but in 
the ideal case itself, where no such deviations are present, 
the existence of ion cores is enough to effectively nullify 
the Peierls instability. 

III. SUMMARY 

To summarize, we have first explored the standard 
half-filled one-dimensional free electron system from the 
perspective of a Peierls distortion and have shown that 
the divergence in the real part of the susceptibility, 
caused by Fermi surface nesting, is exceedingly fragile. 
Effects such as temperature, imperfect nesting, or scat- 
tering, all of which are expected to be operative in a 
real material, reduce the divergence to a simple peak, 
often not more than a factor of two or three enhanced 
over the baseline susceptibility at q~0. Thus, expecta- 
tions of a CDW transition driven entirely by a Fermi 
surface nesting in any real material are unrealistic from 
the outset. Next, we have examined a system in which a 
gap at the Fermi surface has already been opened due a 
commensurate or incommensurate lattice transition, and 
found that the energy gain comes largely from the lower- 
ing of already filled states located away from the Fermi 
energy and not from removing states from the Fermi en- 
ergy itself. This reinforces our contention that the Fermi 



surface topology plays at best a secondary role in CDW 
formation. Using first principles calculations, we take 
three examples of well-known CDW materials to illus- 
trate our point and in each case we find that Fermi sur- 
face nesting either does not exist at the CDW wave vec- 
tor, or is not the strongest nesting in the system. We 
further find that the CDW instability is not fundamen- 
tally electronic, but rather stems from strong electron- 
phonon interaction which, of course, is itself affected by 
the electronic structure.. Finally, we examine the canoni- 
cal Peierls system, a one-dimensional chain of atoms, and 
find that the expected dimerization along the chain axis 
with any realistic amplitude is energetically unfavorable. 
The expected doubling of the unit cell occurs only if the 
one-dimensionality is relaxed to allow a zigzag, rather 
than dimcrized, distortion. We conclude that no true 
distinction between CDWs and structural phase transi- 
tions, in particular incommensurate lattice transitions, 
can be made. Inspecting the Fermi surface itself for possi- 
ble nesting features, without actually calculating the real 
part of electronic susceptibility has no predictive power 
for such structural transitions. Calculating the real part 
of the electronic susceptibility may be helpful in analyz- 
ing such instabilities, but only in relatively few cases can 
it be accepted as the only or even the main driving force 
for such transitions. 
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